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ll i ABSTRACT 

O ; 

^3 ', We present the results from the OGLE-II survey (1996-2000) towards the Large 

Magellanic Cloud (LMC), which has the aim of detecting the microlensing phenomena 
11 ' caused by dark matter compact objects in the Galactic Halo (Machos). 

We use high resolution HST images of the OGLE fields and derive the correction 
y—{ • for the number of monitored stars in each field. This also yield blending distributions 

t> ' which we use in 'catalogue level' Monte Carlo simulations of the microlensing events 

in order to calculate the detection efficiency of the events. 

We detect two candidates for microlensing events in the All Stars Sample, which 
translates into an optical depth of 0.43 ± 0.33 x 10~ 7 . If both events were due to 
Macho the fraction of mass of compact dark matter objects in the Galactic halo would 
' be 8 ± 6 per cent. This optical depth, however, along with the characteristics of the 

events, seems to be consistent with the self-lensing scenario, i.e., self-lensing alone 
. is sufficient to explain the observed microlensing signal. Our results indicate a non- 

detection of Machos lensing towards the LMC with an upper limit on their abundance 
in the Galactic halo of 19 per cent for M — 0.4M Q and 10 per cent for masses between 
0.01 and 0.2 M . 

rS ' 

i— i ' Key words: Cosmology: Dark Matter, Gravitational Lensing, Galaxy: Structure, 

Halo, Galaxies: Magellanic Clouds 
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1 INTRODUCTION 

Although we are constantly learning new things about 
our Galaxy, some major questions about its structure 
still remain unanswered. Among them is the composition 
of the Galaxy halo. Solving this mystery has been the 
main goal of the large ph otometric s urvey s launched in 
the 1990s. The proposal of iPaczvriskil ll 19861) has been re - 
alised b y the surveys like MACHO (lAlcock et alj Il993h . 
OGLE (Udalski et al.lll993l h EROS (|Aubourg et al.1 Il993l ) 
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and MOA (|Yocld ll998l E These groups continuously moni- 
tored millions of stars in both Magellanic Clouds in order 
to detect the transient amplification of brightness caused 
by objects crossing the observer-source line of sight in a 
phen omenon called gravitational microlensing (|Paczvriskil 
1996). The rate of such events depends on the number, mass 
and kinematics of objects populating the line of sight, in the 
disk and the halo of our Galaxy and in the Clouds. The main 
hypothesis all the surveys have been testing is whether the 
halo of our Galaxy contains dark matter in the form of the 
MAssive Compact Halo Objects, called Machos. 

If the Galactic halo were entirely made of Machos, with 
masses between 10 -s and 10 s M@, in a project like OGLE- 
II we should detect about five microlensing events per year 
towards the Large Magellanic Cloud (LMC), and the cor- 
responding optical depth should be TMacho ~ 4.7 x ICE 7 



2 L. Wyrzykowski et al. 



i|Bennettll2005h . However, the actual optical depth measure- 
ments reported by the microlensing experiments are signifi- 
cantly smaller. 

The MACHO collaboration presented a sample of 10 
microlensing event candidates towards the LMC and con- 
cluded that the fraction of dark matter in Machos / = 
MMacho/Mhaio ~ 20 per cent, suggesting the existence of a 
new, previ ously unknown population of objects with masses 
~ 0.5M© l|Alcock et al.l l2000bl . iBennetd 120051 ). Their esti- 
mate of the optical depth towards the LMC is tlmc = 
(1.0 ± 0.3) x 10 -7 . However, the EROS group reported a 
null result and derived an upper limit on the optical depth 
due to M achos of TMacho < 0-36 x 10 -7 , which implies / < 8 
per cent (jTisserand et al 1l2007t ). MACHO's sample of events 
is likely to have some contamination with va riable stars, 
novae, and supernovae ijBelokurov et al.ll2003h . and addi- 
tional data indeed shows the presence of small contamina- 
tion (jBennett et al.l [2005h . However, the optical depth we 
quote f rom MACHO has been corrected for this contami- 
nation |Bennettll2005h . The large discrepancy between the 
results may be due to the different samples of stars used by 
these groups or the fact that they probe different regions of 
the LMC. However, it is unclear whether these fully explain 
the discrepancy. 

Independent searches for Machos have also been car- 
ried out towards the Andromed a Galaxy (M31) usin g pixel - 
lensing method. For example, ICalchi Novati et alj (|2005l ) 
found 6 events and estimated the Machos halo mass frac- 
tion to be high e r tha n 20 per cent. On the other hand, 
iBelokurov et all i|2005h presented an automatic procedure 
that selects unambiguous microlensing events in a pixel light 
curve data set with a total of only 3 "gold-plated" candi- 
dates, thus emphasizing the lack of large nu mbers of clear- 
cut m icrolensing events towards the M31. Ide Jong et alj 
(2006) showed 14 events towards M31 and noticed that the 
observed event rate is consistent with the self-lensing pre- 
dictions, however their results were still inconclusive, mainly 
due to low statistics. 

In this paper we present the results of the search for mi- 
crolensing events towards the Large Magellanic Cloud in the 
independent data set gathered by the Optical Gravitational 
Lensing Experiment (OGLE) in its second phase during the 
years from 1996 to 2000. The paper is organized as follows. 
First, the observational data are described along with the 
error bar correction procedure. Next, the event search al- 
gorithm and its results are presented. It is followed by the 
description of the analysis of blending and its impact on 
the optical depth measurement and the event's detection 
efficiency. Finally, the optical depth is calculated and the 
results are discussed. 



2 OBSERVATIONAL DATA 

The data used in this paper were collected during the sec- 
ond phase of the OGLE survey, in the years 1996 - 2000 
towards the Large Magellanic Cloud (LMC). The OGLE 
project used its own dedicated Warsaw Telescope in the Las 
Campanas Observatory, Chile. D etails of the in s trume nta- 
tion for OGLE-II can be found in lUdalski et all (|l997l ). 

In brief, a 2kx2k CCD with a pixel size of 0.417 arcsec 
was operated in the drift-scan mode, giving an actual size for 



each observed field of 2kx8k. There were 21 fields observed 
towards the LMC covering a total of 4.72 square degrees. 
Their locations are shown in Fig. [T] Fields are listed in Table 
[1] which contains the coordinates of centers of the fields, 
the number of good objects in the /—band, the blending- 
corrected number of stars (see Section and the blending 
density group. By 'good', we mean all objects having at least 
80 epochs during the entire time span of the OGLE-II (from 
about 15 to 25 per cent of all collected frames of a field) 
and mean magnitude brighter than 20.4 mag. The limiting 
magnitude was chosen to be at the peak of the observed 
luminosity function. 

The images in the I and V bands were de-biased and 
flat-field corrected "on-the-fly" . The photometric pipeline 
was based on the Difference Image An alysis method 
fDIA. IWoznTakl |2000| ; TAlard fc Luptonlll998l '). Photometric 
databases were cr eated for both pass-bands as described in 
ISzvmariskil <|2005h and are available on- linfl 

The OGLE-II observations of the LMC began in De- 
cember 1996 (HJD=2450446), but during the first year 
only the very central fields (i.e. those that are, in general, 
denser) were observed (fields 1-10 and 12). In November 
1997 (HJD=2450726), fields 11 and 13-20 were added to the 
observing queue, and the addition of the field LMC_SC21 in 
January 1998 (HJD=2450831) eventually formed the entire 
set of OGLE-II fields. All 21 fields were monitored continu- 
ously until November 2000 (HJD=2451874), yielding about 
500 and 300 frames per field in the / band for dense and 
sparse fields, respectively. In the V band, there were no more 
than 50 frames per field collected in total during the du- 
ration of the project. On average, each field was observed 
every third night in the / band, and every 11-th night in 
the V band. Mean seeing was 1.35 arcsec in the I and 1.37 
arcsec in V. Template images were created by stacking the 
best seeing images, resulting in images with seeing of about 
1.1 arcsec. The template images contained about 5.5 million 
objects suitable for our study in the / and V band. 

The detection of microlensing events and the efficiency 
determination have been performed using I— band data 
only, as these were far more numerous and were sampled 
more frequently as compared to V— band photometry. 

In this study we also occasionally made use of OGLE- 
III LMC data from years 2001-2008, which covered a much 
wider area and included the OGLE-II fields. Details of the 
instrum e ntatio n and photometric pipeline can be found in 
lUdalskil (|2003l ). Analysis of OGLE-III LMC data will be 
presented separately in a forthcoming paper (Wyrzykowski 
et al., in preparation). 

Additionally, for the two detected events, we also use 
the data gathered by the MACHO group and provided on 
their WWW interfaoQ The reductions of the images and 
extraction of the photometry were performed with Differ- 
ence Image Analysis following iKozIowski et all l|2007h . 



3 ERROR-BARS CORRECTION 

Photometric error bars calculated using the Difference Im- 
age Analysis method are known to be generally underes- 

1 http://ogledb.astrouw.edu.pl 

2 http:/ /wwwmacho.mcmaster.ca/Data/MachoData.html 
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Figure 1. Positions of the OGLE— II LMC fields (red rectangles with labels). Also shown are all OGLE-III (green squares) and MACHO 
project fields (blue dashed squares). EROS fields are not shown here as they cover nearly the whole picture. The two small gold-filled 
squares show the positions of the HST fields used for our blending determination. Background image credit: Peter Ward. 



Table 1. The OGLE-II LMC fields. 



Field 






No of stars 


Stellar density 


Density level 








template 


estimated real 


[stars / sq. arcmin] 




LMC.SC1 


5:33:48.84 


-70:06:07.9 


288029 


668799 


367 


dense 


LMC.SC2 


5:31:17.17 


-69:51:57.0 


332267 


771301 


424 


dense 


LMC.SC3 


5:28:47.95 


-69:48:05.6 


355075 


822913 


453 


dense 


LMC.SC4 


5:26:17.82 


-69:48:06.1 


390484 


904432 


498 


dense 


LMC.SC5 


5:23:47.79 


-69:41:03.3 


379709 


878472 


484 


dense 


LMC.SC6 


5:21:18.34 


-69:37:10.3 


392781 


909775 


501 


dense 


LMC.SC7 


5:18:48.01 


-69:24:09.9 


388918 


900339 


496 


dense 


LMC.SC8 


5:16:17.91 


-69:19:13.7 


349202 


807890 


445 


dense 


LMC.SC9 


5:13:47.92 


-69:14:05.8 


289391 


671459 


369 


dense 


LMC.SC10 


5:11:16.10 


-69:09:17.3 


258706 


599375 


330 


dense 


LMC.SC11 


5:08:41.06 


-69:10:03.1 


253395 


588432 


323 


dense 


LMC.SC12 


5:06:16.15 


-69:38:18.3 


186951 


318893 


238 


sparse 


LMC.SC13 


5:06:14.18 


-68:43:30.3 


214229 


365917 


273 


sparse 


LMC.SC14 


5:03:48.62 


-69:04:44.3 


222569 


380394 


284 


sparse 


LMC.SC15 


5:01:17.27 


-69:04:42.9 


170657 


291859 


218 


sparse 


LMC.SC16 


5:36:17.62 


-70:09:41.9 


247981 


575953 


316 


dense 


LMC.SC17 


5:38:47.47 


-70:16:44.1 


194661 


332187 


248 


sparse 


LMC.SC18 


5:41:17.77 


-70:24:48.0 


166998 


285126 


213 


sparse 


LMC.SC19 


5:43:47.91 


-70:34:42.7 


156353 


267126 


199 


sparse 


LMC.SC20 


5:46:17.85 


-70:44:50.0 


152186 


260515 


194 


sparse 


LMC.SC21 


5:21:13.73 


-70:33:18.7 


143979 


244667 


184 


sparse 


total 






5534521 


11845824 







Coordinates point to the centre of the field, each being 14' X 56'. The number of 'good' objects in the template is provided (number of 
data points > 80 and mean /-band magnitude < 20.4 mag) together with the estimated number of real monitored stars (see Section [BJ. 
Stellar density in number of stars per square arc minute was used to classify fields into dense or sparse classes with the threshold of 300 
stars /sq. arcmin. 
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Table 2. Error correction coefficients for OGLE-II LMC fields. 



not corrected 



Field 


7 


e 


LMCJ3C1 


1.327934 


0.002309 


LMC.SC2 


1.534320 


0.000000 


LMC.SC3 


1.514598 


0.000000 


LMC.SC4 


1.534963 


0.000000 


LMCJ3C5 


1.552337 


0.000000 


LMC.SC6 


1.512235 


0.000000 


LMCJ3C7 


1.546070 


0.000000 


LMC.SC8 


1.376516 


0.001565 


LMC.SC9 


1.432857 


0.001066 


LMC.SC10 


1.306555 


0.001539 


LMC.SC11 


1.084232 


0.002695 


LMC.SC12 


1.325315 


0.001657 


LMC.SC13 


1.369874 


0.001256 


LMC.SC14 


0.932450 


0.002931 


LMC.SC15 


1.168966 


0.002046 


LMC.SC16 


1.151916 


0.002628 


LMC.SC17 


1.149929 


0.002365 


LMC.SC18 


1.205946 


0.002624 


LMC.SC19 


1.128117 


0.003138 


LMC.SC20 


1.203747 


0.002623 


LMC.SC21 


1.301942 


0.002467 



timated. This effect is most pronounced for bright stars, 
as noted already by the authors of the DIA method 
(|Alard fc Luptonll 19981 ). In microlensing studies this affects 
the brightest parts of events, increasing the x values in mi- 
crolensing model fitting. In order to increase the detectabil- 
ity of events, we derived and applied an error bar correction 
factor. 

Typically, error bars in microlensing surveys are cor- 
rected by calculating one single correction factor for an 
entire light curve, derived at the event's baseline (as- 
suming the baseline is constant). This approach, however, 
does not account for the need for varying correction fac- 
tors for different apparent magnitudes, which in case of 
highly magnified events leaves the error bars at the peak 
of the event significantly un derestimated. Also, as shown 
in IWvrzvkowski et al] ((2006), some events exhibit intrinsic 
variability in their baselines due to variability of the mi- 
crolensed source or one of the blended stars, in which case 
the derived photometric error bars are overestimated. 

In this study, we processed the entire OGLE-II data 
set and determined the correction coefficients in an em- 
pirical way. For each OGLE-II field, we extracted the I- 
band light curves of stars and compared their intrinsic error- 
weighted rms with the mean error returned by the photom- 
etry pipeline. In order to produce statistically comparable 
results, in each light curve we selected only the first 250 mea- 
surements for further analysis, as different fields had differ- 
ent numbers of frames taken. In the plot of rms vs magnitude 
(see Fig. [SJ), the lowest boundary of the rms is occupied by 
the least variable (i.e., constant) stars. On the other hand, 
the lowest outline of mean error bar is defined by the errors 
assigned to the best measurements at given magnitude. In 
an ideal situation these two curves should follow each other, 
i.e., every constant star should have its rms equal to its mean 
error bar. As can be seen from the upper panel of Fig. [21 for 
our data the mean errors were systematically lower than the 
rms. 
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Figure 2. Mean error (green) and rms (black) for a sample of 
stars from OGLE-II LMC-SC1 field before (upper panel) and 
after (lower panel) the correction of error bars. Lines indicate the 
lower outlines of the rms (red) and mean error bar (blue). 



In order to account for this systematic shift we have 
fitted the following two-parameter formula which minimised 
the difference between lower outlines of the rms and mean 
error bars: 



AI cor = V(7A/) 2 + e 2 



(1) 



where Al is the error bar returned by the photometry 
pipeline, 7 and e are the correction coefficients and AI cor 
is the corrected error bar. 

The lower panel of Figf2] shows the rms and mean er- 
ror bars after applying the corrections to each light curve's 
individual error bar measurements using Equation [T] 

Table [5] presents 7 and e coefficients derived for each 
OGLE-II LMC field. The mean 7 and e for all LMC fields 
were 1.32 and 0.0016, respectively. Derived correction coef- 
ficients are suitable for use in any study involving OGLE-II 
Large Magellanic Cloud photometric data. 

As a by-product of the error bar analysis we derived 
the formula for scaling real error bars from some reference 
light curve to any magnitude, which is necessary for the 
simulation of light curves: 

■,0.35(J slm -J ref ) 



AJ sin 



A7 ref 10 



for J sim > 15.0 



(2) 
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Figure 3. Two sample light-curves for false-positive events, 
which were rejected in the search procedure. Top: "blue bumper" 
(which exhibited a second bump in the OGLLE-III data). Bot- 
tom: possible nova or supernova with additional variability in the 
declining part. Solid magenta line is the best microlensing model, 
dashed blue line is a nova model. 



0.003, for I uim < 15.0, 



where J s i m is the simulated magnitude for which the error 
bar (A/sim) is required, I le f and A7 re f are the magnitude 
and the error bar of the reference star at a given epoch. 
These error bars still need to be corrected using eq. ([l]). 



4 SEARCH PROCEDURE 

We searched for microlensing events among 5.5 million ob- 
jects, selected from the database as having more than 80 
data points and mean magnitude in I band brighter than 
20.4 mag (cut 0). This defined our All Stars Sample. All 
light curves now have their measurements' error bars cor- 
rected as described in Section [3] 

Table [3] presents all of the cuts, along with the number 
of objects left after each cut. All parameters for this proce- 
dure were derived and fine-tuned for OGLE-II LMC data 
after carrying-out Monte Carlo simulations of microlensing 
events, using the simulation procedures described in Section 

ra 

The first step of the search procedure is similar to the 
one used in the search f or microlensing ev ents in the OGLE- 
II Galactic bulge data (|Sumi et al.ll2006l ). All objects were 



searched for a bump over a baseline, allowing for some vari- 
ability in the baseline. 

For a given light curve, we introduce an outer window 
with width equal to half of the total observational time-span. 
The significance of each data point Oi with respect to the 
outer window was calculated via 

Oi = , = (3) 

where Ii is the magnitude of the point i and Ali is its error 
bar, while I m ed,w and a w are the median and the rms in the 
outer window, respectively. The bump was detected when 
there were at least 4 consecutive data points above a a > 1.4 
threshold and 3 of them where above a a > 2.5 threshold. 
Up to 2 peaks were allowed in one light curve, but then only 
the one with the highest a max was chosen. If the sum of 
as over the peak data points was higher than 30.0 the light 
curve was selected (cut 1). 

In cut 2. we m asked out possible "blue bumpers" 
l|Alcock et al. 2000b) - variable stars located on the bright 
blue part of the Main Sequence, exhibiting outbursts that 
can look similar to the microlensing phenomenon. The up- 
per panel of Fig. [3] shows an example of a "blue bumper" 
star, rejected due to its magnitude and colour. This partic- 
ular object also exhibited additional significant variation in 
the OGLE-III data in subsequent years. 

In the next step, the rema i ning l ight curves were fit- 
ted with the standard IPaczvhskil (|l996l ) microlensing model, 
(i.e., a point-source - point-lens microlensing event) which is 
described as: 



/ = Jo -2.5 log [/ s A+(l-/s)] : 
where 



A = 



u z + 2 
tiVu 2 + 4 



and 



(4) 



(5) 



The fitted parameters are: io - the time of the max- 
imum of the peak, £e - the Einstein radius crossing time 
(event's time-scale), tto - the event's impact parameter, 7o - 
the baseline magnitude in the I band and fs - the blending 
fraction (ratio of lensed source flux to total blends' flux in 
the I band). The fits were performed in two ways, namely 
with blending parameter fixed fs = 1, i.e., with no blend- 
ing, and with fs being free. For clarity, parameters of the 
non-blended (4 parameters) model are given the subscript 
M 4. 

Then we compared the 4-parameter microlensing model 
(with blending fixed to 1 to assure convergence) with a con- 
stant line fit (cu t 3). A similar cut wa s also applied in the 
EROS analysis l|Tisserand et al.ll2007l ). where it proved to 
be very efficient. Ndof,n4 and Ndof,iU,peak are the degrees of 
freedom of the fit calculated using numbers of data points in 
the entire light curve and around the peak, respectively. This 
cut removed many low signal-to-noise bumps and allowed us 
to select only the outstanding bumps. 

In cut 4, we required that there were more than 6 data 
points around the microlensing peak in the range of io±te/j4. 
This meant that we had enough data points in the bump to 
analyse. 

Next we attempted to remove all plausible novae, su- 
pernovae and all kinds of asymmetric bumps contaminating 
our sample (cut 5). We fitted the light curves with an gen- 
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Table 3. Selection criteria for microlensing events in the OGLE-II data 



Cut no. 



No. of objects left 



Selection of "good" objects 

1 Significant bump over baseline 

2 "Bumper" cut^ 

3 Microlensing fit better than constant line fit 

4 Number of points at the peak* 

5 Microlensing fit better than supernova fit, 

6 Peak within the data span 
[HJD-2450000] 

7 Blended fit converged 

8 Conditions on goodness of microlensing fit 
(global and at the peak) 

9 Time-scale cut 

H 

10 Impact parameter cut 



N > 80, (7) ^ 20.4 mag 

V" <Tj > 30.0 

peak 

(I) > 19.0 mag, (V - I) > 0.5 mag 
x , 2 x »"'~*£ 4 > 140 

«do/ 4 M 4 V 2N dof.„4,peak 

X 2 SN >MIN( x \xld 
446 < t ^ 1874 



< fs < 1.2 



if— < 2.3 and x ^-^ k <; 5 



1 < t E <C 500 



< M ^ 1 



5,534,521 
9,390 

7,867 
1,020 

968 
607 
591 

113 
4 
2 



t magnitudes as in the field LMCJ3C1 (shifted according to the position of the center of Red Clump) 

*in the range of to ± tE^4 



eral asymmetric mo del, following t he fo rmula used for the 
novae search in e.g. iFeenev et al.l l|2005t ). which comprises 
of two scalable exponents and a constant baseline. We com- 
pared its goodness-of-fit with the x 2 of microlensing model, 
taking the smaller of blended and non-blended fits. This cut 
was sensitive to all sorts of asymmetric light curves. Among 
the removed objects were several dwarf novae candidates, 
as well as long-term varying objects and several eruptive Be 
stars. These Be stars were fainter or redder than our Blue 
Bumper cut and typically e xhibited a flat bas eline, a sharp 
rise and a slow decline (e.g. jKeller et al1l2002l ). 

Assumin g a supernova rate of 0.5 SNe yr _1 deg -2 above 
V = 20 mag l|Alcock et all lioOOb), then for our experiment 
of 4 years covering 4.5 square degrees and assuming detec- 
tion efficiency of about 20 per cent we should expect about 
1 supernova in our data set. Visual inspection of the light 
curves rejected at the supernova fit stage (cut 5) indeed 
showed one plausible candidate for a supernova (see lower 
plot of Figure [3]l . 

In the following step (cut 6) , we removed all light curves 
in which the peak occurred before the beginning or after the 
end of the time-span of the OGLE-II observations. 

Next, in cut 7 we selected only those candidates whose 
blended microlensing fit has converged with blending pa- 
rameter fs < 1.2. We allowed for small am ount of 'negative ' 
blending followin g previous studies (e.g. JPark et all |2004 , 
ISmith et al.ll2007l ) showing regular event can be affected by 
this effect. The vast majority of the light curves removed in 



cut 7 had fs significantly larger that 2, indicating that the 
blending fit did not converge. 

In cut 8, we requested the reduced x 2 of the blended 
microlensing fit was less than 2.3 and the \ 2 °f the non- 
blended model calculated on the peak only (within 1 tE^4 
around to) was less than 5.0. This stage removed obvious 
variable stars with a single bump lying above the rest of the 
light curve to which the microlensing model was fit quite 
poorly. We were able to apply this cut uniformly to all light 
curves, independently on their magnitude, thanks to the fact 
that the error bars were corrected and any dependence on 
the measurement brightness was accounted for (Section [3j. 

The final cuts 9 and 10 are applied to select candidates 
suitable for the optical depth determination. For this, we 
used microlensing parameters obtained with the fitting pro- 
cedure, and required the time-scales tE to be between 1 and 
500 days and impact parameter uo less than 1 Einstein ra- 
dius. 



5 RESULTS 

Out of 5.5 x 10 6 light curves we found only two plausible can- 
didates for microlensing events. These are labeled OGLE- 
LMC-01 and OGLE-LMC-02. 

Table [4] presents their coordinates, OGLE-II field, star 
id, baseline I— and l/-band magnitudes obtained in the mi- 
crolensing fit and magnitude and colour of the source. 

Both candidate light curves were checked for avail- 



The OGLE View of Microlensing towards the Magellanic Clouds. I. OGLE-II LMC 7 



Table 4. Microlensing candidates detected in the OGLE— II LMC data. 



Event's name 


RA 


Dec 


field 


database 


baseline / 


baseline V 


source I 


source V — I 




[J2000.0] 


[J2000.0] 




star id 


[mag] 


[mag] 


[mag] 


[mag] 


OGLE-LMC-01 


5:16:53.26 


-69:16:30.1 


LMC.SC8 


235928 


19.91 


20.65 


19.90 


0.887* 


(EWS: OGLE-1999-LMC-01) 










±0.01 


±0.14 


±0.06 


±0.007 


(MACHO-99-LMC-2) 


















OGLE-LMC-02 


5:30:48.00 


-69:54:33.6 


LMC.SC2 


170259 


20.42 


20.68 


20.53 


0.46 












±0.02 


±0.13 


±0.38 


±0.03 



* colour derived using MACHO B-band data 




Figure 4. Light curves of the two OGLE-II events OGLE-LMC-01 (left panel) and OGLE-LMC-02 (right panel). Black points are 
the OGLE / band, blue triangles are the MACHO B-band and OGLE V-band observations, for OGLE-LMC-01 and OGLE-LMC-02, 
respectively. Solid red and green lines show the best-fitting microlensing models to the data, with residuals shown below each curve. The 
inset shows a zoom-in around the peak of amplification, indicating that the finite-source effect is not detectable. Re-reduced MACHO 
data and the fit for OGLE-LMC-02 are shown above the OGLE data shifted by 0.5 mag (squares for B, crosses for R). 



able additional data outside the time-span of the OGLE- 
II project. We cross-correlated them with the OGLE-III 
database, covering years 2001-2008 (the most-recent data 
used were obtained in May 2008) and the MACHO database, 
covering years 1992-1999, yielding total data coverage of 16 
years. Neither candidate exhibits additional variability apart 
from the actual event. 

In the OGLE-II database there were 44 l/-band observa- 
tions of OGLE-LMC-01, but none of them were taken during 
the event. Therefore for colour determination of the source, 
we used re-reduced MACHO data in the B band. By com- 
paring the OGLE V and I CMD with MACHO B-band pho- 
tometry of nearby stars, we found the following transforma- 
tion for the colour: (-B-J) = (1.1571±0.0046)(V-/)±const, 
where the measured constant is an artefact of the instru- 
mental magnitude system of our reductions of the MACHO 
data. For OGLE-LMC-02, we used the OGLE V-band data, 
as there were 59 observations available, with about 14 taken 
during the peak. 

Close-ups of the events are shown in Fig. U along with 
the best-fitting microlensing models and their residuals. 

The top panel of Fig. [5] shows the colour-magnitude 
diagram of LMC st ars based on OGLE and HST 
l|Holtzman et al.l |2006J ) data, along with the determined 
positions of the sources in the two OGLE candidates. 



For comparison, we have also included in this figure 
the candidates from the MACHO collaboration: differ- 
ent colours and shapes differentiate b e tween self-lensing 
candidates according to iMancini et al.l l|2004l ). the binary 
event MACHO-L MC-9, the thick disk lens candidates 
MACHO- LMC-5 feallivavalil et all 12004 ) and MACHO- 
LMC-20 (|Kallivavalil et al l 12006) and the remaining can- 
didates. 

In the bottom panel of Fig. [5j the positions of OGLE- 
II and MACHO candidates are shown on the map of stellar 
density of the Red Clump stars from the OGLE-III data. 



5.1 OGLE-LMC-01 

This was the only event detecte d towards the LM C by 
the Early Warning System (EWS, llJdalski et al.|[l99ll ) dur- 
ing the entire second phase of the OGLE project and 
was designated OGLE-1999-LMC-OI0. It was also inde- 
pendently discovered and alerted by the MACHO collab- 
oration (MA C HO-99 -LMC-02), but it was not included in 
lAlcock et al.l l|2000bl ) as it occurred outside the 5.7 year 



3 http: / / www.astrouw.edu.pl / ~ogle/ogle2 / ews /1999 /lmc- 
01. html 
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Table 5. Parameters of the microlensing model fits to the OGLE-LMC-01 and OGLE-LMC-02 events. 



parameter 


OGLE-LMC-01 
5-parameter fit 4-parameter fit 


7-parameter fit 


to 


1337.20 


±0.02 


1337.20 


±0.02 


1337.20 


±0.01 


*E 


57.2 


+4.7 
-4.2 


59.2 


±0.7 


65.7 


±2.6 


«o 


0.0258 


+0.0024 
-0.0023 


0.0248 ±0.0004 


0.0231 


±0.0009 


Io 


19.91 


±0.01 


19.91 


±0.01 


19.92 


±0.01 


fSi 


1.04 


±0.09 


1.0 




0.91 


±0.04 


Bo 










18.28 


±0.01 












1.11 


±0.05 


2 

x 


286.7 




286.9 




3401.0 




x 2 

N dof ■■■■ 


0.88 




0.87 




2.68 




parameter 


OGLE-LMC-02 
5-parameter fit 4-parameter fit 


7-parameter fit 


to 


491.6 


±0.3 


491.6 


±0.3 


491.3 


±0.3 


tE 


23.8 


+6.0 
-5.1 


24.2 


±1.0 


27.4 


+6.9 
-5.6 


UO 


0.4120 


+0.1856 
-0.1164 


0.4011 


+0.0083 
-0.0081 


0.3478 


+0.1428 
-0.0955 


Io 


20.42 


±0.02 


20.42 


±0.02 


20.44 


±0.02 


fSz 


1.04 


+0.77 
-0.37 


1.0 




0.83 


+0.52 
-0.28 


V 










20.68 


±0.02 


/s v 










0.67 


+0.42 
-0.22 


x 2 


578.2 




578.2 




636.2 




x 2 

N daf ■■■■ 


1.21 




1.20 




1.19 





time-sp an of analysed dat a (MACHO star ID: 79.5863.4522, 
see also iBond et al.1 12002| ) . 

The OGLE-II data were fitted with the standard 
Paczyhski microlensing model, i.e., a point-source-point-lens 
model (eq. [4} . When we fixed blending such that the event 
was unblended (/sm4 = 1) we obtained £e/j4 = 59.2 ± 0.7 
days, whilst when fs was free, we found t E = 57.2± 4 4 7 2 days 
(fs = 1.04 ± 0.09). When OGLE I- and MACHO B-band 
data were fitted simultaneously, the derived timescale was 
somewhat larger (£e vi = 65.7 days), as the blending param- 
eter in I-band changed slightly. The results of the fits are 
summarised in Table [5] 

Since the event was relatively long we checked 
if O GLE-LMC-01 was subject to the parallax effect 
fe.g. jHan fc GouTd| [l995b'), The components of the parallax 
vector 7te in the OGLE-II data fit were 7te,jv = — 0.03±0.10 
and 7te,b = —0.16 ± 0. 25, i.e., consis tent at 1-cr level with 
no detection of parallax (Gould 2004). Therefore, the time- 
scale derived in the fit is not subject to any degeneracy with 



parallax. In principle the lack of parallax signatures allows 
us to place limits on lens distance. However, for this event 
the parallax constraints are not very tight and so the limits 
are not particularly insightful: if we assume that the source 
lies at the distance of the LMC (~ 50kpc) and the lens is 
O.5M , we find that the 2-sigma lower-limit on the lens dis- 
tance is 300 pc. From this, we cannot make any significant 
statements about the population to which the lens belongs. 

With a maximum amplification of about 40 (~ 4 mag), 
this is one of the most secure and best constrained events to- 
wards the LMC found to date. High amplification events are 
often prone to exhibit finite-source effects in cases where the 
source and the lens are close to each other. We checked for 
the presence of this effect in the OGLE and MACHO data. 
In particular, MACHO data had two observations very close 
to the peak, where the effect should be most pronounced. We 
have not detected any deviation (see inset in Fig. [4}. This 
only allows us to put an upper limit on the source star size 
in units of the Einstein radii of p > 0.0332, which along 
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with an estimation of the source size {9» = 0.42 ^tas) lead to 
an upper limit on — 9 t /p > 13.1 /jas. This at the LMC 
distance is equivalent to about 0.66 AU. Hence, the limit 
on the projected transverse velocity of the lensing object 
v± = 0E/tE > 21 km/s, which is satisfied by typical lenses 
in both the LMC and Galaxy halo. 

A simultaneous microlensing model fit to the OGLE I- 
band and MACHO B-band data gives an apparent source 
magnitude of I s = 19.90±0.06 and colour {V-I) a = 0.887± 
0.007. Magnitudes and location of the source on the CMD 
in Fig. [5]indicate it belongs to the Red Giants Branch of the 
LMC. 



5.2 OGLE-LMC-02 

The second event, OGLE-LMC-02, was found just at the 
limit of our magnitude cut (the actual mean brightness was 
20.29 mag, which means that it passed cut 0). The maxi- 
mum amplification of ~2.6 produced a rise in brightness of 
~1 mag. This event occurred in the very first season of the 
LMC observations in the OGLE-II, when the sampling was 
much denser than in subsequent seasons. Therefore, the mi- 
crolensing light curve was very well covered and prominent 
(see Fig. 0}. 

Again, we used the standard point-source-point-lens 
model to fit the OGLE-II data. The unblended fit (/ S(l4 = 1) 
gave a time-scale £e^4 = 24.2 ± 1.0 days, and for uncon- 
strained f s we found t E = 23.8^1 days (/ s = 1M±%]£). 
The large error bars in fs and tv were also present in the si- 
multaneous fit to I- and V-band observations, which yielded 
a time scale of t-E lv ~ 27 .4 days. 

This object was also observed by the MACHO collabo- 
ration (MACHO star ID: 77.8152.2917). However, the event 
was not detected by them, probably due to noisy data and 
sparse sampling. The peak was clearly revealed in both red 
and blue filters only after DIA re-reductions of the original 
MACHO images and its best-fitting microlensing model had 
a time-scale of around 25 days, consistent with the fit to the 
OGLE data. 

From the fit to the two-band OGLE data, we derived 
the magnitude and colour of the lensed source. The position 
of the source with magnitude Is = 20.53 ± 0.38 and colour 
(V — I)s = 0.46±0.03, is marked on the Fig.[5]and indicates 
it is probably a Main Sequence star located in the LMC. 

Its location on the Main Sequence branch of the CMD 
means that this may potentially be a candidate for a 
blue bumper. However, while many blue bumpers exhibit 
secondary brightening episodes, the baseline of this event 
proved to be constant over more than 15 years. Therefore, 
we believe that it is unlikely to be a blue bumper. 



6 BLENDING AND NUMBER OF 
MONITORED STARS 

Microlensing collaborations have deliberately chosen to 
monitor the densest fields on the sky, namely the Galac- 
tic Center and Magellanic Clouds. This is done in order to 
increase the number of potentially detectable microlensing 
events, as the chances of seeing an events are of the order 
of 10~ 7 (i.e. , one out of 10 million monitored stars is under- 
going a microlensing episode at any given time). However, 




-0.5 



0.0 



0.5 1.0 
V-l [mag] 



1 .5 



2.0 




RA [2000] 

Figure 5. CMD of OGLE-II and MACHO candidate mi- 
crolensing events, overploted on the OGLE (red background) and 
HST (black dots) measurements of stars in central parts of the 
LMC (upper panel). Positions of OGLE and MACHO candidate 
events on the OGLE— III LMC Red Clump star count map (lower 
pane l). The contours sho w modelled self-lensing optical depth 
from lMancini et al.1 (120041 ') and run every 0.4 X 10~ 8 with 4 X 10~ 8 
at the innermost one. Dashed boxes mark OGLE-II fields. Posi- 
tions of the sources in the two OGLE events are marked in red 
and with a number. Remaining symbols represent MACHO can- 
didates with binary event #9 (green star), candidates for self- 
lensing (yellow squares), confirmed thick-disk lenses #5 and #20 
(dark blue filled squares) and remaining marked with blue trian- 
gles. 



in such crowded stellar fields observed with a medium-sized 
ground-based telescope, virtually every object detected on 
the image consists of several stars blended together. Neglect- 
ing blending can have se vere consequences on the optical 
depth measurements fe.g JSumi et al.| [2006). as the ampli- 
fications and time-scales of microlensing events may be de- 
termined incorrectl y, sometimes by a factor of 2 or more. 
On the other hand, ISmith et al.l |2007) showed that in the 
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case of the Galactic bulge the effect of overestimation of the 
events' time-scales cancels out with the underestimate of the 
real number of monitored stars to lowest order. 

Dealing with blending is a delicate and difficult task. 
Adding a blending parameter to the standard microlensing 
light curve model may cause degener acies in the parameter 
space of the microlensing events fe.g. JWozniak fc Paczvriskil 
119971 ). Therefore, the favoured approach in recent optical 
depth determinations towards the Galactic bulge and Mag- 
ellanic Clouds is to neglect blending by narrowing the star 
sample to the bright est stars only |Popowski et alj 120051 . 
iTisserand et alj|2007ft . However, in the case of Magellanic 
Clouds the number of possible events is very small and it is 
even more unlikely to detect a microlensing event occurring 
in the limited number of bright stars. 

In our analysis, we decided to use all stars available in 
the OGLE-II database. To deal with the blending issue - 
which is most serious for faint stars - we make use of the 
archival images from the Hubble Space Telescope (HSTj3 
located by chance in our fields. The HST images have supe- 
rior resolution when compared to the ground-based OGLE 
images (i.e. they are nearly blending-free) and are usually 
much deeper than the OGLE data. Therefore these data are 
a valueable source of information regarding the amount of 
blending. In the HST archive, we identified two deep ACS 
/-band images coinciding with the OGLE fields, located in 
OGLE fields LMC.SC6 and LMC.SC19, as typical repre- 
sentatives of the dense and sparse parts of the LMC. The 
location of these fields is shown in Fig. [T] 

The positions of the HST stars were transformed onto 
the OGLE images and then we identified which HST stars 
contributed to each OGLE object by calculating their dis- 
tances from the OGLE blend's centre. A blend's radius was 
estimated using the Gaussian profile for the given blend's 
magnitude at a chosen level above the background (7 ~ 19.7 
mag) with an additional 1.3 pixel added in quadrature. We 
identified HST stars that were not more than 3.5 mag fainter 
than the OGLE blend (i.e.,/s > 4 per cent) and consid- 
ered stars located within the blend radius as components of 
that blend. This method allowed us to create an empirical 
distribution of the blending fractions for dense and sparse 
fields. We divided these distributions into three magnitude 
bins: 14 mag< I ^17.5 mag, 17.5 mag< I ^19 mag and 
19 mag< I ^20.4 mag. They are shown in Fig. [6] 

Then, every OGLE LMC field was classified as dense or 
sparse according to its average stellar density down to 20.4 
mag (with 300 stars per square arcmin as a division bound- 
ary) on the OGLE template image (see Table [1} • We then 
adopted the corresponding distributions in the detection ef- 
ficiency calculations and in our recalculation of the number 
of monitored stars (see Section [7J. 

Due to blending in ground-based observations, we 
can not resolve individual stars that are separated by less 
than the typical seeing (of order of 1 arcsec). Even though 
the OGLE template images, with respect to the typical im- 
ages, have usually high resolution thanks to stacking of ~a 
dozen of the best-seeing images, we still have limited abil- 
ity to resolve individual stars. Since each OGLE template 
object may consist of several stars, any of which could be 



4 http://archive.stsci.edu/hst/ 
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Figure 6. Distributions of blending fractions for stars from the 
analysis of the archival high-resolution HST images of parts of 
the OGLE fields. 
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number of HST stars in one OGLE star 

Figure 7. Distribution of the number of HST stars in one OGLE- 
II object obtained for dense and sparse fields. 



microlensed, the number of stars we consider as potential 
sources for microlensing events is actually larger than num- 
ber of objects detected on the template image. 

In order to estimate the actual number of monitored 
stars - which enters directly into eq. (J6]l - we derived the 
distributions of the number of HST stars hidden in a sin- 
gle OGLE template object (see Fig. using the same HST 
images and photometry. Then, for each OGLE field we con- 
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volved its luminosity function with the distribution in each 
magnitude bin (for a corresponding density level). The esti- 
mated numbers of monitored stars in each field are provided 
in Table [T] On average, the correction factor for dense fields 
was about 2.3 and for sparse about 1.7. The total number of 
monitored stars was estimated to be about 11.8 million, com- 
pared to about 5.5 million objects detected on the OGLE-II 
template images in cut 0. 

It is worth pointing out that in our method of count- 
ing the total number of monitored stars we did not count all 
stars that would be observed with the HST down to 23.9 mag 
(20.4 mag + 3.5 mag). We counted only those stars which 
were blended with other stars such that they formed a blend 
of magnitude brighter than 20.4 mag on the OGLE template 
image. For instance, a single 21 mag star might not be de- 
tected on the template image. Even if microlensed, we would 
not detect this event as in our analysis we did not search 
for transient events with no baseline (i.e., with no object on 
the template). If such star, however, was blended with an- 
other 21 mag star they would form a 20.2 mag object - very 
likely detectable on the template. In this case, both compo- 
nents would have blending factor of /s = 0.5 and both of 
them would be counted as a separate potential source for 
microlensing events. 
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Figure 8. Detection efficiencies for OGLE-II LMC.SC8 field for 
All (thin lines) and Bright (thick lines) Stars Samples, with blend- 
ing included (solid lines) and neglected (dashed lines). For most 
of the dense fields, the detection efficiencies follow these relations. 
The efficiencies for the most of the sparse fields are lower by a 
about 30 per cent due to the fact they have a shorter time-span 
of observations. 



7 DETECTION EFFICIENCY 

The task of calculating the detection efficiency for microlens- 
ing events can be a complicated procedure. In general, this 
requires us to simulate events and then determine the frac- 
tion that we recover with our selection criteria. It may seem 
that the most accurate method is to inject microlensing 
events into the series of real frames or to create completely 
new mock images and then put them through the entire pho- 
tometric pipeline. However, this is obviously very intensive 
computationally. Furthermore, it is not clear how to deal 
with the fact that any star within the blended object can be 
microlensed or whether the mock frame is indeed reproduc- 
ing all characteristics of the real one. 

Some of these problems can be also encountered when 
simulating microlensing events at the so-called 'catalogue- 
level' (which only involves generating series of light curves). 
However, the computational effort required for this proce- 
dure is significantly reduced and some difficulties can be 
resolved more easily. 

The 'image-level' simulations were performed by MA- 
CHO and MOA collaborations in their calculation of the op- 
tical depth towards the Galactic bulg e (jAlcock et al. 2000a; 
ISumi et al.ll2003l ; lAlcock et al.ll200lh . The 'catalogue-level' 
appr oach has been more commonly used for the bulge 
(e.g.jAfonso et alJliooll ; lHamadache et al.ll2006l ; ISumi et all 
l200fj ) and in the EROS collabor ation's study of microle nsing 
towards the Magellanic Clouds (|Tisserand et al.ll2007l ). 

In this study, we decided to use the existing photomet- 
ric light curves as a base for the simulations. However, as 
we discuss below, we have attempted to incorporate all of 
the advantages o f the 'image-level' simulations. As shown in 
iKozIowskil l|2007h such an approach can be very successful 
when all factors are carefully taken into account. 

The detection efficiency was determined individually for 
the fields containing both events (SC2 and SC8). For a wide 



range of time-scales between 1 and 1000 days, we performed 
Monte Carlo simulation of about 10 events. From the field's 
photometric database, a star was picked randomly, provided 
that the star satisfied the basic cut (i.e., it had the required 
number of data points and sufficient mean brightness). For 
a given £e, the other microlensing parameters were drawn 
from flat distributions (to and uo) and from the empirically 
derived distributions for the blending parameter (/s; see Sec- 
tion |5J). The latter was chosen with respect to the stellar 
density level of the simulated field and the magnitude bin 
in which the simulated star was located. The original flux 
of the input star was apportioned between the source and 
blend in the ratio fs : 1 — /s- 

Effectively, the lensed flux was added to already existing 
stars in the database, so no new stars were 'created'. Also, 
the procedure preserved any variability and non-gaussianity 
present in the original photometry, allowing for effects of de- 
creasing or increasing of the variability amplitude during mi- 
crolensing (|Wvrzvkowski et al.l I2006T ) . Photometric errors- 
bars were scaled according to the empirically derived rela- 
tion between the error bars and the magnitude (see eq. [3}. 
Such simulated events were then analysed with our search 
pipeline described in Section [4] 

In their LMC data, the MACHO collaboration found 
one event (MACHO-LMC-9) to be a clear caustic crossing 
binary lens. Our detection pipeline discriminated against 
binary and other exotic events, as we only fitted a point- 
source-point-lens microlensing model. However, we per- 
formed a visual inspection of the several thousands of high 
signal-to-noise light curves that passed cut 1, and found no 
candidates for any non-standard type of events. Our detec- 
tion efficiency did not take into account the influence of 
binary lenses because Monte Carlo simulations of binary 
events and their automated recovery is a difficult and com- 
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Figure 9. Upper panel: observed luminosity functions in two 
OGLE-II LMC fields, where the dotted and solid lines represent 
sparse (SC19) and dense (SC6) fields, respectively. The vertical 
line shows the cut-off at I = 20.4 mag. Lower panel: luminosity 
functions for the same fields recovered after applying the blend- 
ing correction (see Section [6}. Their shapes follow the prototype 
LFs of two HST fields (dense and sparse) used for the blending 
correction (thick dash-dotted lines). 



putationally demanding task. Theref ore we followed the ap - 
proach adopted by the EROS group l|Tisserand et al.|[2007l ) 
and reduced the efficiency by a factor of 0.9 to compensate 
for the 10 per cent possible binary lenses. 

We performed our simulations on both Bright and All 
Stars Samples with blending included. In addition, we re- 
peated the analysis assuming no blending (i.e., fs — 1) 
in order to compare our results with other studies. We de- 
fined the Bright Stars Sample as stars with N > 80 and 
brighter than 7rc + 1 mag, where 7rc is the Red Clump 
stars' centre on the colour-magnitude diagram. In the refer- 
ence field LMC.SC1, this threshold was at 18.8 mag. The All 
Stars Sample comprised of all stars used in the search for 
microlensing events, i.e., stars brighter than 20.4 mag and 
with at least 80 data points. 

Fig. |8]shows the derived detection efficiencies for a wide 
range of time-scales for a typical dense field. Most sparse 
fields were observed for a shorter period of time (obser- 
vations started about a year later than for dense fields). 
However, in our detection efficiency calculations we simu- 
lated microlensing events within the full time-span of the 
OGLE-II experiment, namely for to within the range of 



HJD = 2450446 - 2451874. Therefore, the efficiency for 
fields observed for a shorter time is reduced by about 30 per 
cent, owing to the events which fall outside the simulated 
light curve. Only one field classified as sparse (LMC_SC12) 
was monitored since the very beginning of the project and 
the microlensing detection efficiency for this field is compa- 
rable to the efficiencies for dense fields. 

Fig. [8] also shows the detection efficiencies for the non- 
blended analysis of All and Bright Stars Samples (with fs 
1 for all simulated events). The efficiency for Bright Stars 
with and without blending is higher than the corresponding 
efficiencies for All Stars Sample, as it is generally easier to 
recover events occurring on brighter stars. 

The observed luminosity functions (LFs) of dense and 
sparse OGLE-II LMC fields are shown in the upper panel of 
Fig. Blending is more severe in dense fields and it causes 
faint stars to 'merge' and form brighter objects, populat- 
ing brighter parts of the luminosity functions. The bottom 
panel of Fig. [5] shows blending- free luminosity functions for 
a dense and sparse field, recovered using the blending distri- 
butions obtained by comparing OGLE and high-resolution 
HST images. The luminosity functions obtained from HST 
are shown as well. It is noticeable that the effect of over- 
populating the brighter parts and under-populating the faint 
end in blending-free LFs of dense fields has virtually disap- 
peared after applying the correction. The recovered dense 
and sparse luminosity functions follow their prototypes from 
the HST everywhere except at the bright and very faint end. 
For stars brighter than I ^ 17 mag the HST images be- 
come saturated, which is reflected in the steeper luminosity 
functions at the bright end. The large scatter in the bright 
range is due to low number statistics in the HST's small 
field-of-view (as compared to OGLE-II fields). The broad- 
ening of the Red Clump may also be caused by this effect, or 
could be due to the distance gradient along the OGLE field. 
The HST's LFs of dense and sparse areas differ significantly, 
which is likely caused by the different mix of LMC bar and 
disk populations in these two fields. 

In our Monte Carlo simulations, we drew an object from 
the original OGLE luminosity functions. Therefore we were 
not affected by the incompleteness of the bright end of the 
HST luminosity functions. Incompleteness in the faint end 
of the OGLE luminosity functions is partially recovered in 
the correction process, but also is negligible as the sensitiv- 
ity for such faint stars is very low. Although our blending 
analysis successfully recovers the underlying stellar bright- 
ness distribution, it is still approximate as we used only two 
HST images to represent dense and sparse fields. In future, 
it would be interesting to expand this study with a detailed 
7— band high-resolution follow-up, using HST and/or other 
high-resolution instruments. 



8 OPTICAL DEPTH ESTIMATE 

We determined the optical depth using the All Stars Sample. 
For the two events found in the OGLE-II data we used 
the following standard equation for calculating the optical 
depth: 

iVov 
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Table 6. The optical depth for the two OGLE-II events. 



event 



n x icr 7 



efficiency not corrected for binary events 



j+4.7 



OGLE-LMC-02 23.8^1° 0.165536 0.13 



total 7XMC-02 



0.38 ±0.29 



efficiency corrected for binary events 



OGLE-LMC-01 57.2+l- 7 2 0.190957 0.28 



-,+6.0 



0.15 



total r LMC _ 2 



0.43 + 0.33 



where T b s = 1428 days is the time-span of all observations, 
N* = 11.8 X 10 6 is the total number of monitored stars (see 
Section [6]), -/V cv = 2 is the total number of events, tEi is the 
time-scale of each event detected with the efficiency of e(tEi)- 
For the time-scales of these events we used the microlens- 
ing fit where blending was a free parameter (5-parameters 
model). The error in r was c alculated with the formula 
given bv lHanfc Gouldl | 1995al ). For the original detection 
efficiency (i.e., not corrected for binary events), the optical 
depth was derived to be tlmc-02 = 0.38±0.29xl0~ 7 . If the 
efficiency is corrected for non-detectability of binary lenses, 
the optical depth becomes tlmc-02 = 0.43 ± 0.33 x 10 -7 . 
Table [6] presents all calculations of the optical depth. 



9 DISCUSSION 
9.1 Optical depth 

The OGLE-II LMC data analysed in this study covered 
about 4.5 sq. deg. of the sky over 4 years, which is clearly 
smaller than the data used by MACHO or EROS groups. 
However, it provides a completely independent data set to 
test the hypotheses concerning the presence of MACHOs in 
the Galactic halo. 

Based on two the events we found the optical depth 
towards the LMC to be t L mc-02 = 0.43 ± 0.33 x 10" 7 . 
On one hand, this value is in 2-a agr eement with t he 
Tlmc = 1-0 ± 0.3 obtained by MACHO l|Bennettl I2005M 
If both OGLE-II events were caused by Machos, then they 
would contribute 8±6 per cent to the total mass of the halo, 
according to model S from lAlcock et al.l l|2000bl ) . 

On the other hand, our value of the optical depth is 
more consistent with a scenario of lensing due to the visi- 
ble stellar LMC component alone (self-lensing) , for which 
various estimates of the optical de pth are present in the lit- 
erature: r SL = (0.05-0.78) x 10~ 7 (|Gvuk et al.ll2000T l. r SL = 



5 Only 5 (#1, #8, #14, #15 and #18) out of 17 MACHO can- 
didates occurred within OGLE— II fields. However, all of these 
events peaked before the first observations of OGLE— II. There- 
fore we can only confirm continuity of their constant baselines 
using our data. 



0.24-0.36) x 1(T 7 jAlcock et al.ll200o"bh , r SL 
Belokurov et al.ll2004h.~ 



0.54 x 10" 



Mancini et ail (|2004h presented a detailed theoretical 



study of the halo and self-lensing optical depth and com- 
pared it with the results of MACHO. They simulated the 
disk and bar of the LMC and quantitatively derived a con- 
tribution from the LMC's own stars to the optical depth as 
r < 6 x 10~ 8 . In the region of the LMC bar, where most 
OGLE-II fields are concentrated, they estimated the self- 
lensing optical depth to vary between 2 and 4 x 10~ 8 . The 
optical depth derived for two events we found in the OGLE- 
II data and averaged over all fields is in a perfect agreement 
with these estimates. 

Also the time-scales of our two events seem to belong 
to the regime of self-lensing, especially the 57.2-days event, 
whic h lies exactly on the li ne for the self-lensing events in fig. 
10 of lMancini et alj <|2004h . From the CMD of the LMC, it is 
apparent that both microlensed sources are coincident with 
an LMC population, indicating that they probably belong 
to the LMC. The locations of both events in the bar of the 
LMC also supports the self-lensing scenario. 

If we assume that our two events are indeed due to self- 
lensing (i.e., we have zero events caused by the halo objects), 
we can derive an estimate for the upper limit of the Macho 
abundance in the halo of the Milky Way. We performed this 
analysis f ollowing the pro c edure similar to the one used by 
EROS in lTisserand et all l|2007l ). first estimating the num- 
ber of expected events in the OGLE-II data s et if the halo 
was en tirely made of Machos, using model S |Alcock et al.l 
|2000bJ). In this model the optical depth towards the LMC 



equals 



= 5.1 x 10" 



which also includes a contribu- 



tion from self-lensing and Galaxy disk lensing. Dark mat- 
ter lensing solely had Tjjjacho = 

4.7 x 10 -7 . For a so-called 
'typical' mass of 0.4 Mq, in OGLE-II we should expect 16 
events. At the 95 per cent Poisson confidence level there 
are up to 3 events consistent with a zero detection. This di- 
vided by a number of expected events gives the upper limit 
on the mass fraction of about 19 per cent. This translates 
to an upper limit on the optical depth due to Machos of 
toam q < 0.9 x 10" 7 . If blending is neglected in the All Stars 
Sample the numbers are somewhat larger: about 9 expected 
events would put a limit of ro.4Af Q < 1.6 x 10~ 7 . 

Fig. (TTJshows the dependence of the number of expected 
events on the mass of the lensing objects for the blended 
and non-blended approach. Fig. Illl shows the 95 and 90 per 
cent confidence limits on the halo Macho mass fraction for a 
wide range of masses. It also compares our results with the 
results from EROS and MACHO. This figure clearly shows 
that the OGLE-II project is the most sensitive to events 
caused by Machos with masses between 0.01 and 0.2 Mq. 
For this range of masses, we can place an even tighter upper- 
limit of about 10 per cent (at 95 per cent confidence) on the 
abundance of compact dark matter objects in the Galactic 
halo. A more detailed theoretical analysis will be presented 
in a forthcoming paper (Calchi Novati, Mancini, Scarpetta 
et al., in preparation). 

For the sake of comparison with the results of the EROS 
group, we also performed the above analysis narrowing the 
stars sample only to the bright stars. In each field, our 
Bright Sample consisted of about 20 - 30 per cent of stars 
present in All Stars Sample, with a total over all fields of 
about 1.9 million. After correction for blending, the num- 
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logM = 2log(<t E >/70d) logM = 2log(<t E >/70d) 



Figure 10. The number of microlensing events caused by Machos 
expected to be seen in the OGLE— II data, for blended and non- 
blended approaches. 



ber of monitored stars was about 3.6 million. In the Bright 
Sample of OGLE-II, we have n ot found any events . Ac- 
cording to the full S model of lAlcock et ail (|2000bh . we 
should expect about 9 events to be caused by lenses of 
mass 0.4Mq, including events due to self-lensing. If we ne- 
glect blending, then we expect 6 events. The total opti- 
cal depth towards the LMC is therefore estimated to be 
toaMq < 1.6 x 10~ 7 (r .4Af Q < 2.6 x 10~ 7 ) with (without) 
the full blending correction. The calculated upper-limit is 
very high due to the fact that OGLE-II has much smaller 
sky coverage and shorter time-baseline compared to the com- 
bined EROS-I and EROS-II projects. However, a combined 
study of OGLE-II and OGLE-III data will be conducted in 
the near future and should be able to yield constraints of 
even strong er magnitude than those given by the EROS col- 
laboration l|Wvrzvkowski et alj|200sl Wyrzykowski et al., in 
preparation). 

9.2 Blending 

As has been discussed above, we did not have HST data for 
all OGLE-II fields. Therefore we chose to use the two se- 
lected HST fields and consider these as being representative 
of two different density levels. In this way we recover only 
an approximation to the underlying real luminosity function 
for each field. 

In order to check how our selection of HST stars af- 
fects the result, we repeated the whole analysis for a dif- 
ferent cut on the HST stars magnitude, namely 2.5 mag 
(instead of 3.5 mag) below the brightness of the OGLE ob- 
jects. This translates to a minimal allowed blending fraction 
of 10 per cent. We then rederived the blending distributions 
and the distributions for the number of HST stars in a single 
OGLE object. These distributions were obviously different 
compared to the original 3.5 mag cut; for both sparse and 
dense fields, there was nearly no blending in the brightest 
magnitude bin. For the fainter stars, the distributions indi- 
cated a level of blending similar to that which was found 
previously, but with slightly less pronounced features for 
small fs- 

In a similar manner to that described in Section [6] we 



Figure 11. Mass fraction in compact dark halo objects as a 
function of the mass of the lensing objects. The red box with a 
star shows the value for the two OGLE-II events assuming they 
are caused by halo lenses. Solid and dashed lines show OGLE-II 
upper limits assuming both OGLE events are due to self— lensing. 
Also shown are the results of the MACHO collaboration (dotted 
line) and the upper limit derived by the EROS group (dot-dashed 
line). 



estimated the total number of monitored stars and found 
about 9.2 x 10 6 stars. Then we also repeated the efficiency 
calculation including these new blending distributions. It 
was similar to the original efficiency to within 2 per cent. 
The detection efficiencies (corrected for binary events) for 
the two events were found to be e(tE = 57.2) = 0.214 and 
e(tE = 23.8) = 0.165. Since these efficiencies are similar to 
previous values while the number of monitored stars is re- 
duced, we find that this yields a somewhat larger optical 
depth of THST2.5 = (0.49 ± 0.37) x 10~ 7 . This is, however, 
still consistent with our original value, indicating little sen- 
sitivity of the optical depth on the depth of blending stars' 
magnitude limit. 

As an extreme case we also considered neglecting blend- 
ing completely in our All Stars Sample. In this case, we find 
that there are only about 5.5 million sources with fs = 1. 
The detection efficiency was obtained as described above 
and is shown on Fig. [5] It was larger than in the case when 
blending is taken into account by about 20 per cent. How- 
ever, this does not cancel out entirely with the reduction in 
the number of monitored sources. Using these numbers with 
non-blended model time-scales (Table [5} leads to a larger 
value for the optical depth of r = (0.72 ± 0.55) x 10~ 7 , 
where we have used the detection efficiencies corrected for 
lack of binary events. 

This indicates that neglecting blending favours larger 
values of the optical depth, emphasising the importance of 
carrying out a careful analysis of blending when dealing with 
All Stars samples in order to obtain reliable optical depth 
results towards the LMC. 



10 CONCLUSIONS 

In the search for microlensing events in OGLE-II data to- 
wards the Large Magellanic Cloud, we found two events. 
The total optical depth for all 21 OGLE-II LMC fields, cov- 
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ering about 4.5 square degrees, was derived for the All Stars 
Sample (with a limit of (J) < 20.4 mag) as tlmc-02 = 
0.43±0.33 1(T 7 . If both events were caused by Machos, this 
would imply that their comprise 8 ± 6 per cent of the halo 
by mass. 

However, the value of the optical depth and the char- 
acteristics of the events are consistent with the self-lensing 
scenario, in which lensing towards the LMC occurs only be- 
tween sources and lenses both located in the Cloud. From 
this, we can constrain the presence of compact dark matter 
objects in the Milky Way's halo and we find that our results 
are compatible with no detection. The upper limit is found 
to be 19 per cent for Muacho = 0.4Mq and 10 per cent for 
masses between 0.01 and 0.2 M@. It is worth emphasising 
that in our study we probed mainly the bar of the LMC, 
similarly to the MACHO project. However, our results and 
conclusions are closer to the ones of the EROS collaboration, 
which probed different regions of the LMC. 

The sensitivity of the OGLE-II experiment is limited 
to a relatively short time span and sky coverage. However, 
the results presented in this paper will be soon strengthened 
and supplemented with the analysis of the 8-year OGLE-III 
data. 

Future photometric and astrometric surveys may con- 
firm the lensing origin of most of events detected so far and 
the new ones by high-resolution imaging or precise astrome- 
try, allowing for measuring the distances to the sources and 
lenses or their velocities. 
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